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Abstract. We present a 100 yr optical lightcurve 
of the recently discovered supersoft X-ray source 
RX J0019. 8+2156 as deduced from photographic plates of 
Harvard and Sonneberg Observatory. Apart from the pe- 
riodic orbital variations two different timescales of optical 
variability are discovered. The timescales and amplitudes 
of this variability are discussed in the framework of the 
steady nuclear burning model. 
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1. Introduction 

After the discovery observations of supersoft X-ray sources 
(SSS) with Einstein , the ROSAT satellite has discovered 
two dozens new SSS (Greiner et al. 1991, Schaeidt et al. 
1993, Greiner et al. 1994, Motch et al. 1994, Kahabka 
et al. 1994, Beuermann et al. 1995, Supper et al. 1995). 
The ROSAT spectra of SSS are well described by black- 
body emission of very low temperature (kT« 25-40 eV) 
and a luminosity close to the Eddington limit (Greiner et 
al. 1991). Several SSS exhibit X-ray variability. The first 
of the new ROSAT discovered SSS - RX J0527.8-6954 
- must have been a factor of ten fainter during the ear- 
lier Einstein observations when it was in the FOV but 
not detected. A transient SSS was observed during the 
ROSAT All-Sky-Survey to rise from below 0.06 cts/sec to 
about 2 cts/sec within a few days (Schaeidt et al. 1993). 
Later observations showed the source being alternately on 
and off on a few months timescale (Schaeidt et al. 1995). 
These observations suggest that the variability occurs on 
timescales which are short compared to the nuclear burn- 
ing timescale of novae or the limit cycle length as proposed 
by van den Heuvel et al. (1992). 

RX J0019. 8+2156 was found during the search for 
galactic SSS in the ROSAT All-Sky-Survey data and op- 
tically identified with a V=12.2 mag object (Beuermann 



et al. 1995). Its optical spectrum is very similar to that of 
CAL 83 (Cowley et al. 1984). As the optically brightest 
SSS, RX J0019. 8+2156 is best suited for a study of the 
long-term optical behaviour of this class of objects. 

2. Observational results 

2.1. Photographic plates 

The Harvard plate collection comprises plates taken at dif- 
ferent observatories and with different instrumentations 
(catalogued as a specific series each). Plates of the series 
I, RH, A, AC, MC, and Damon were used for this study 
(see Table 1). From the Sonneberg plate collection the blue 
sensitive plates of the sky patrol (SHU) have been used. 
They were gained in general with Zeiss Ernostar and Tes- 
sar type cameras of focal lengths between approximately 
20 cm and 30 cm. The majority of plates used here are 
centered at /l +20°. Our comparison stars (Fig. 1) have 
been linked to the system of the Mt. Wilson International 
magnitudes of Selected Area No. 68 (Seares et al. 1925). 
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Fig. 1. Finding chart of RX J0019. 8+2156 (cross in the center) 
with the photographic magnitudes of the comparison stars. 
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Fig. 2. Optical lightcurve of RX J0019. 8+2156 derived from photographic plates at Harvard (lozenges) and Sonneberg (triangles) 
Observatory, respectively. Each data point represents the mean of one season. Filled symbols contain more than 9 individual 
measurements with la statistical errors being less than 0.05 mag, whereas open symbols contain 9 measurements or less. The 
dashed box marks the part shown in more detail in Fig. 3. The lower panels are blow-ups of the two time intervals covering the 
sudden intensity jumps. Here, the data are shown in bins of three months together with la error bars. 



The use of different emulsions in the past 100 years 
together with the existence of strong emission lines in 
the optical spectrum of RX J0019. 8+2156 demanded 
some care in the data reduction. Thus, plates with un- 
known emulsion type or filter combination have been dis- 
carded. In addition, several other stars in the vicinity 
of RX J0019. 8+2156 have been monitored to ensure the 
recognition of rare, falsely documented information on 
emulsion and/or filter. We estimate the error of a single 
magnitude determination to be 0.15 mag. 

2.2. Optical variability 

In addition to the optical variability already described in 
Beuermann et al. (1995), namely 

1. a cyclic variation of a period of 15.8 hours with an 
amplitude of about 0.3 mag, 

2. quasiperiodic "pulsations" of roughly 2 hours length 
each and a range of <0.1 mag, 



we detected two further components of variability: 

3. long-term variations, seemingly non-periodic, with 
timescales of up to 20 years and an amplitude of 1 
mag (Fig. 2) and 

4. irregular fluctuations of timescales of weeks to months 
and small range (Fig. 3). 
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Fig. 3. Typical variation on the timescale of weeks. The large 
symbol represents 9 magnitude determinations, the remaining 
are single measurements (all Sonneberg plates). Data of the 
phase interval 0.9-0.1 have been discarded. Arrows mark the 
mean of the preceding and following seasons. 
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Table 1. Archival plates used 



Plate series m.u m (mag) Time span No. of plates 
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Fig. 4. Orbital variation of RX J0019. 8+2156 in three time 
intervals. The period (as given in the text) is stable at least 
over the last 40 years. 

There are two large optical brightness jumps by about 
1 magnitude. The rise time of the jump in 1969 is shorter 
than 10 months, while the jump in 1931/32 has a rise time 
of 14 months. Both of these intensity jumps are followed 
by a slow relaxation continuing over «20 years. Whether 
or not there was a similar jump at the end of the last 
century (which would imply a nearly periodic recurrence 
time scale of «40 years) cannot be decided due to the 
sparse coverage. But certainly, at these times the object 
also showed the full amplitude of variability (>1 mag). 

Many of the seasonal means as plotted in Fig. 2 are av- 
erages over irregular fluctuations as shown in Fig. 3. These 
variations occur on timescales of weeks to months, and the 
variability pattern are different from year to year. Thus, 
a considerable part of the scattering in Fig. 2 can be as- 
signed to these irregular, short-term fluctuations with the 
resulting mean being dependent on the different coverage 
of bright and faint episodes. 

In the best-covered sections of the Sonneberg material, 
when also emulsion types and cameras used were homo- 
geneous, component 1 can be easily traced despite of the 
other components. Our observations present no contradic- 
tions to the period given by Beuermann et al. (1995). We 
were, however, able to substantially improve it and derived 
the following lightcurve elements (ephemeris): 

Min.(heL) = 243 5799.247 + 0'?6604565 x E 
which is valid at least for the years 1955-1993 (Fig. 4). 
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3. Discussion 









3.1. Binary dimension 

The stability of the period over many years is most easily 
explained by orbital motion. The very broad minimum, 
and the obvious vanishing of the "pulsational" feature 
(above item 2) during the central parts of the minimum 
point towards an eclipse of a bright accretion disk. The 
presence of a strong reflection effect which would show the 
same phase behaviour as an eclipse phenomenon, cannot 
be excluded by our observations alone, but is not proba- 
ble in view of the spectral behaviour (Beuermann et al. 
1995). The dimension of the circular orbit - half axis a = 
P 2 / 3 (Mi+M 2 ) 1/3 -fits well with the half amplitude of the 
observed radial velocity curve - v = 27r((Mi+M2)/P) 1 ' 3 
giving a « 4 solar radii, M1+M2 ~ 2 solar masses. We 
conclude therefore that one of the most promising mod- 
els for the object is an accreting white dwarf (WD) with 
nuclear burning in a classic cataclysmic system. We note 
that models comprising a high mass component or a giant 
secondary star (classic symbiotic binary) can be excluded. 

3.2. The long-term lightcurve 

To variable star experts long-term lightcurves with 
timescales of several years and amplitudes of the order of 
1 mag are not uncommon. Shell stars of the Pleione type 
(see Hoffmeister et al. 1985 for composite lightcurves) or 
slowly varying red giants and supergiants (ibid) exhibit 
such behaviour and the objects of the present kind could 
be mistaken for them if only photometric observations 
were at hand. However, we also know various long-term 
phenomena in mass transferring binaries. Examples are 
long standstills between eruption phases (Z Cam) or in 
outburst (UX UMa), inactive states in LMXBs (HZ Her), 
long lasting bright or faint episodes in "nova- like" vari- 
ables (KR Aur, MV Lyr, TT Ari), or a slowly varying 
minimum brightness in dwarf novae (V 426 Oph). In most 
of these cases a long-term change of the mass transfer 
emerging from the secondary is the simplest explanation. 
RX J0019. 8+2156 had an X-ray intensity of of 2.1 
cts/sec during the All-Sky-Survey (1990), and 2.0 cts/sec 
in pointed ROSAT observations in 1992 and 1993 (Beuer- 
mann et al. 1995), i.e. there was no variation in X-ray 
luminosity over 3 years. Also, no temperature changes 
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could be found. During the same time interval the mean 
optical intensity has increased by nearly 0.5 mag. Unfor- 
tunately, there are no exactly simultaneous ROSAT and 
optical measurements. Given the short-term variability on 
timescale of weeks (Fig. 3), we cannot draw a firm conclu- 
sion on the relation of X-ray and optical luminosity. 

3.3. Burning white dwarfs 

The most popular model involves steady nuclear burning 
on the surface of an accreting WD (van den Heuvel et al. 
1992). Hydrogen burning on the WD surface is stable for 
M > 10" 7 M©/yr (Paczynski & Zytkow (1978). However, 
the burning timescale (100-300 yrs) is much longer than 
the observed optical variability of RX J0019. 8+2156. Pos- 
sibly, the observations can be understood by temporar- 
ily reduced mass transfer which leads to long-term "os- 
cillations" between stable and unstable H burning. Con- 
sequently, the optical brightness jumps by 1 magnitude 
could be caused be ignitions of H (or He) burning: Ex- 
tensive steady state calculations of massive white dwarfs 
accreting at high rates have shown (Iben 1982, Fujimoto 
1982) that for WD masses in excess of 1.0 M@ the on/off 
timescales of H flashes become shorter than 50 years for 
accretion rates of the order of 10 -7 M@/yi. Hydrodynamic 
calculations have shown (Livio et al. 1989) that the lumi- 
nosity in the off state increases with increasing M, thus 
resulting in smaller outburst amplitudes. Thus, one could 
imagine that (possibly varying) mass transfer near the 
burning rate drives the WD only into mild flashes. Af- 
ter having burned the rather small hydrogen layer, the 
system returns to instability. In order to avoid large am- 
plitude flashes, the subsequent cooling (along the constant 
radius track in the HR diagram) has to be stopped by re- 
suming accretion at a high rate. Thus, hydrogen ignition 
on massive WDs would not necessarily be connected with 
large amplitude intensity variations. 

An interesting complication in the above scenario is 
that though no mass is ejected at accretion rates of about 
10" 7 M©/yr the WD expands to «80 R @ (Livio et al. 
1989) being much larger than the binary dimension. Sion 
& Starrfield (1994) therefore proposed to consider ex- 
tremely hot but low-mass WDs (« 0.5-0.7 M@) which re- 
main at compact radii (R-wd < 0.1 R @ for 10 -8 M @ /yr). 
However, it remains unclear how the intensity variations 
and amplitudes as observed in RX J0019. 8+2156 can be 
explained. Certainly, shell flashes in such systems are 
much too violent, and the deduced cycle lengths are a 
few thousand years (ibid). 

The above scenario has attributed the cause of the 
variability to changing mass transfer rates which in turn 
lead to variable nuclear burning rates and hence depar- 
tures from a limit cycle recurrence time. Consequently, 
one should also ask for the possible reasons of mass trans- 
fer variations. First, if the donor star is more massive than 
the compact object, mass transfer is unstable on a thermal 



timescale. Second, if indeed a F or G type dwarf star (as 
suggested for CAL 83 and CAL 87) is the mass sponsor- 
ing secondary component, our hypothesis is that variable 
median and long-term intrinsic activity of the secondary 
("solar cycle activity", depending on the size of the con- 
vective atmosphere) might give rise to a changing mass 
transfer, and therefore to an optical variability. Spectro- 
scopic investigation of the mass transfer and accretion rate 
during low and high states can clarify this point. Third, 
also intensity and/or spectral variations of the irradiation 
of the secondary can affect the mass transfer. 

4. Conclusion 

In addition to the optical variability of RX J0019. 8+2156 
reported in Beuermann et al. (1995) we have discovered 
two further timescales of variability. We believe these to be 
consequences of unstable mass transfer onto a WD which 
sporadically is burning hydrogen. However, some of the 
details of the proposed models are not in line with the 
observed amplitude and timescale of variability. In addi- 
tion, the applicability of some of the previously applied 
assumptions such as accretion via a disk instead of spher- 
ical accretion, the effect of the expected strong radiation 
driven wind, the angular momentum loss carried by the 
escaping gas and tidal effects in these tight binaries might 
be explored more carefully. 
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